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Abstract 

%^J . A phenomenological parameterization of the proton polarized structure function gi{x, Q'^) 

\^ [ has been developed for x > 0.02 using deep inelastic data up to ~ 50 (GeV/c)'^ as well 

C^ ■ as available experimental results on both photo- and electro-production of proton reso- 

ld . nances. According to the new parameterization the generalized Drell-Hearn-Gerasimov 

(^ I sum rule is predicted to have a zero-crossing point at Q^ = 0.16ib0.04 (GeV/c)'^ . Then, 

low-order polarized Nachtmann moments have been estimated and their Q'^-behavior 
has been investigated in terms of leading and higher twists for Q^ > 1 (GeV/c)'^ . The 
leading twist has been treated at NLO in the strong coupling constant and the effects 
of higher orders of the perturbative series have been estimated using soft-gluon resum- 
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TT*- mation techniques. In case of the first moment higher-twist effects are found to be 

^ . quite small for Q^ > 1 (GeV/c)'^, and the singlet axial charge has been determined to 

be ao[10 (GeV/c)'^] = 0.16 it 0.09. In case of higher order moments, which are sensitive 
to the large-x region, higher-twist effects are significantly reduced by the introduction 
of soft gluon contributions, but they are still relevant at Q^ ~ few (GeV/c)'^ at variance 

5^ I with the case of the unpolarized transverse structure function of the proton. Our find- 

ing suggests that spin-dependent correlations among partons may have more impact 
than spin-independent ones. As a byproduct, it is also shown that the Bloom-Gilman 
local duality is strongly violated in the region of polarized electroproduction of the 
A(1232) resonance. 
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1 Introduction 

The experimental investigation of lepton deep-inelastic scattering (DIS) off proton and 
deuteron targets has provided a wealth of information on parton distributions in the nu- 
cleon, leading to a nice confirmation of the leading (LO) and next-to- leading order (NLO) 
predictions of the perturbative Quantum Chromodynamics (QCD). In the past few years 
some selected issues in the kinematical regions corresponding to large values of the Bjorken 
variable x have attracted a lot of theoretical and phenomenological interest; among them one 
should mention: i) the flavor decomposition of the parton distributions, with particular 
emphasis on the ratio of d- to u-quark as a; — * 1 (see Refs. ||l|, 0), and ii) the occurrence of 
power corrections associated to dynamical higher-twist operators measuring the correlations 
among partons (see Refs. 01-0)- The extraction of the latter is of particular relevance since 
the comparison with theoretical predictions either based on first-principle calculations (like 
lattice QCD simulations) or obtained from models of the nucleon structure may represent 
an important test of the non-perturbative QCD regime. 

Various analyses of power-suppressed terms in the world data on the unpolarized 
nucleon structure functions F^(a;, Q^) and F^(x, Q^) have been carried out in the past 
years. They are based either on the choice of a phenomenological ansatz or on renormalon- 
inspired models [|l], 0, 0, adopting for the leading twist the LO or NLO approximations. 
The effects of the next-to- next-to-leading order (NNLO) corrections have been investigated 
on F.^{x,Q^) and R^,rp{x,Q'^) in Ref. p|, and on the parity-violating structure function 
xF^i^x^Q"^) in Ref. 0. Very recently [§] the effects of high-order radiative corrections on 
the extraction of leading and higher twists in the transverse structure function F^{x, Q"^) 
have been considered adopting soft-gluon resummation techniques 0. It has been shown 
that the extraction of higher twists at large x is remarkably sensitive to soft gluon effects as 
well as to the updated PDG value of as{Mz) 0. Existing analyses indicate that for x > 0.7 
dynamical power corrections in F^{x,Q'^), R^,rp{x,Q'^) and xF^ {x,Q'^) are not very large. 
We want to point out that only in Refs. [|l| and p the Q^ range of the analyses has been 
extended down to Q^ ~ 1 (GeV/c)'^ thanks to the use of Nachtmann moments, including in 
this way the contributions of both the nucleon-resonance regions and the nucleon elastic peak. 
This is clearly worthwhile not only in order to enhance the sensitivity to power-suppressed 
terms, but also because of parton-hadron duality arguments [|| . 

The aim of this paper is to extend the twist analysis made in Refs. [|l| and [|| to 
the case of the polarized proton structure function gi{x,Q'^). As in Ref. |1]], we make 
use of the formalism of the Operator Product Expansion {OPE) [|l^ and of the (polarized) 
Nachtmann moments |]TT|, as defined in Refs. [|r^, |T^, in order to disentangle the kinematical 



target-mass corrections from the dynamical higher-twist effects we are interested in. The 
evaluation of the Nachtmann moments requires however the knowledge of the polarized 
structure functions in the whole x-range for fixed values of Q^. Therefore, we have developed 
a new parameterization of gKxyQ"^), which describes the DIS proton data up to Q^ ~ 
50 [GeV/cY and includes a phenomenological Breit-Wigner ansatz able to reproduce the 
existing electropro duct ion data in the proton-resonance regions. Our interpolation formula 



for gi{x, Q^) has been successfully extended down to the photon point, showing that it 



nicely reproduces the very recent data [|14[ on the energy dependence of the asymmetry 
of the transverse photoproduction cross section as well as the experimental value of the 
proton Drell-Hearn-Gerasimov {DHG) sum rule |jl5|. According to our parameterization 



of gi{x, Q"^) the generalized DHG sum rule is predicted to have a zero-crossing point at 
Q2 = 0.16 ±0.04 {GeV/cf. 

Then, low-order polarized Nachtmann moments have been evaluated and their Q"^ 
behavior has been investigated in terms of leading and higher twists for Q^ > 1 {GeV/cf'. 
The leading twist is treated at NLO in the strong coupling constant and it is extracted 
simultaneously with phenomenological higher-twist terms from our pseudo-data. The effects 
of higher orders of the perturbative series are estimated using the same soft-gluon resum- 
mation technique adopted in case of the analyses of the unpolarized data made in Ref. [0. 
The main results of our power correction analysis are as follows. As far as the first moment 
is concerned, the effects of higher twists are found to be quite small for Q^ > 1 {GeV/cY. 
Moreover, the singlet axial charge is determined to be ao[10 [GeV/cY] = 0.16 ± 0.09; our 
extracted value is significantly below the naive quark-model expectation (i.e., compatible 
with the well-known "proton spin crisis"), but it does not exclude completely a value of the 
singlet axial charge as large as ~ 0.25, in nice agreement with recent estimates (see, e.g., 
Ref. |jl6|). In case of higher order moments, which are more sensitive to the large-z region, 
higher-twist effects are significantly reduced by the introduction of soft gluon contributions, 
but they are still relevant at Q^ ~ few {GeV/c)"^, at variance with the case of the unpolarized 
transverse structure function of the proton (see Ref. 0). Our finding suggests that spin- 
dependent correlations among partons may have more impact than spin-independent ones. 
As a byproduct, it is also shown that the Bloom-Gilman (BG) local duality |]1^ is strongly 



violated in the region of polarized electroproduction of the A(1232) resonance. 

The paper is organized as follows. In the next Section the Nachtmann definition of 
the moments and the NLO approximation for the leading twist are briefly reminded. In 
Section 3 a new parameterization of gi{x,Q'^), which describes the DIS regime as well as 
the photo- and electro-production proton-resonance regions is presented and adopted for the 
evaluation of the Nachtmann moments. Moreover, the issue of a possible local BG duality 
among the DIS behaviour of gi{x, Q"^) and suitable local averages in the resonance regions 
is addressed. Section 4 is devoted to a twist analysis of our pseudo-data at NLO, while the 
inclusion of the effects of high-order radiative corrections is presented in Section 5. Finally, 
our main conclusions are summarized in Section 6. 

2 The Nachtmann moments and the Leading Twist at 

NLO 

The complete Q^ evolution of the structure functions can be obtained using the OPE ||T0[ 
of the time-ordered product of the two currents entering the virtual-photon nucleon forward 



Compton scattering amplitude, viz. 



T[J{z) J(0)] = Y: /n(-^') z>^^z^\..z^- O 
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where O" are symmetric traceless operators of dimension c?" and twist r" = rf" — n, 

with a labehng different operators of spin n. In Eq. (|l]) /°(— z^) are coefficient functions, 
which are calculable in pQCD at short-distance. Since the imaginary part of the forward 
Compton scattering amplitude is simply the hadronic tensor containing the structure func- 
tions measured in DIS experiments, Eq. (P leads to the well-known twist expansion for the 
Cornwall-Norton {CN) moments of the nucleon polarized structure function g^{x,Q'^) (see 
Refs. [0, H), viz. 



T=2,even 



where n = 1, 3, 5, ..., /i is the renormalization scale, Onri^) are the (reduced) matrix elements 
of operators with definite spin n and twist r, containing the information about the non- 
perturbative structure of the target, and EnT{fi',Q'^) are dimensionless coefficient functions, 
which can be expressed perturbatively as a power series of the running coupling constant 

For massless nucleons only operators with spin n contribute to the n-th CN moment 
(0). When the nucleon mass M is taken into account, operators with different spins can 
contribute and consequently the higher-twist terms in the expansion of the CN moment 
M^ {Q"^) contain now also target-mass terms, which are of pure kinematical origin and 



therefore of no physical interest. It has been shown by Nachtmann [|ri[| in the unpolarized 
case and subsequently generalized to the polarized structure functions in Refs. |jl2|, |l3l that 
even when M ^ the moments can be redefined in such a way that only spin-n operators 
contribute in the n-th moment, namely 

Px^ An 
g2 n + 2 



-9f(.,0=)^^l (3) 



where n = 1, 3, 5, ... and 



2x 



1 + 4M2xVQ2 
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is the Nachtmann variable. Note that at variance with the unpolarized case, where the 
transverse Nachtmann moments involve only the transverse structure function F-^lXjQ"^) 
(see, e.g., Ref. ^), in the polarized case MJ^^Q"^) has to be constructed using the two 
polarized structure functions gi{x, Q"^) and 5'^(x, Q"^). Using the experimental data for the 



latter ones in the r.h.s. of Eq. (^, the target-mass corrections are exactly canceled out; 
therefore the twist expansions of the experimental Nachtmann moments MJ^^Q'^) contain 
only dynamical higher twists, which are the only ones related to the correlations among 
partons, viz. 



M^^\Q^) = djj^^^Q"^) + dynamical higher twists 



(5) 



where 5fi^\Q'^) stands for the n-th CN moment of the leading twist contribution. The 
latter can be written in the following form: 
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where at NLO |TB| 
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and 
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In Eq. (^ ^Qn^if^"^) is the ra-th moment of the non-singlet [NS) polarized quark distribution 
in the (modified) minimal subtraction (MS) or Adler-Bardeen (AB) scheme evaluated at the 
renormalization scale /x^, 7^*^ and 7^'^*^ are the non-singlet (unpolarized) anomalous dimen- 
sions at one- and two-loop levels, respectively, SCJ^^ is the n-th moment of the quark coeffi- 
cient function, (3o = ll-2Nf/3, /3i = 102-38iV//3 and eventually < e^ >= (1/iVj) E^S e?, 
with Nf being the number of active flavors. In Eq. (§) (5S„((5^) and SGniQ"^) are the ra-th 
moments of the singlet-quark and gluon polarized distributions, respectively, and 5C^^^ is 
the n-th moment of the gluon coefficient function. 

As it is well known, the evolutions of the singlet-quark and gluon distributions are 
coupled in general; however, as it happens in the unpolarized case (cf., e.g., [Q), the moments 
of order n > 1 are sensitive to the large-x regions, where the evolution of gluons and quarks 
are approximately decoupled. Therefore, in what follows we will assume the following NLO 
evolution for 5/i„>3: 
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where 6An{fi^) = 6fi';^\fi^)/[l + as{fi^)6Ci^y2n]. 



In case of the first moment {n = 1), adopting the AB scheme as defined in Ref. |]T9 
and introducing the simphfied notation A/i^^^ = Sfi^^ii, one gets 



A/i^^HQ 



< e^ > 



Aq^' + aoiQ' 
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where ao(Q^) is the (scale dependent) singlet axial charge, given by 



ao(g') = AS-iV^ 



2tt 



-AG{Q' 



(11) 



In Eqs. (|lOf 11) both the non-singlet Ag^"^ and the quark singlet AS are conserved quantities, 
while the polarized gluon moment AG is scale dependent (at NLO one has AG((5^) ~ 
l/as(Q^), so that the axial singlet charge oq becomes only slightly scale dependent). Below 
the charm threshold (i.e., Nf = 3) one has 



Ag 



NS 
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where qa = 1.2670 ±0.0035 is the nucleon axial coupling constant and ag = 0.579 ±0.025 
^ is the octect axial charge (obtained from nucleon and hyperon beta decays under the 
assumption of S'?7(3)-fiavor symmetry). Thus, the expected value for Aq^^ below the charm 
threshold is Aq^^ = 1.095 ± 0.007. 



3 Construction of the polarized Nachtmann moments 

For the evaluation of the Nachtmann moments M^'{Q^) (Eq. (^) systematic measurements 
of the structure functions gi{x,Q'^) and f7^(a;,Q^) are in principle required in the whole 
x-range at fixed values of Q^. The kinematical coverage of existing data on gi{x, Q"^) is 
shown in Fig. 1. It can be seen that the available data cover the kinematical range 0.3 < 
Q^ (GeV/c)'^ < 60 with values of the produced invariant mass W"^ = M^ ± Q^(l — x)/x up 
to ~ 300 GeV^. While data are scarce below x ~ 0.02, the coverage at intermediate and 
high values of x, but still in the DIS regions (i.e. ior W >^2 GeV), appears to be sufficient 
for developing an interpolation formula. On the contrary, in the nucleon resonance regions 
{W < 2 GeV) only data for Q'^ ~ 0.5 and ~ 1.2 {GeV/c)"^ are presently available from the 
E143 experiment ^1 . 

The polarized proton structure functions 5^1 (x,Q^) and g2{x,Q'^) (we omit the suffix 
p for simplicity) are related to the measured helicity-dependent virtual photon-nucleon cross 
sections by 
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Figure 1. Kinematical coverage of existing data on gi{x,Q'^). Tlie dashed line corresponds to an 
invariant produced mass equal to W^ = 2 GeV. Full dots, open dots, full squares, open squares, full 



diamonds, open and full triangles correspond to the experiments of Refs. |2l|, |2^, |23|, |2^, ^, |26|, |27| 



respectively. Open diamonds are the existing data in the proton-resonance regions from Ref. |21|. 



where ai/2 and o"3/2 are the transverse absorption cross sections for total helicity 1/2 and 
3/2, respectively, and air is the longitudinal-transverse [LT) interference cross section. In 
Eqs. (|T^-|T^) K is the incoming flux factor, which in the Hand convention is explicitly given 
by K = u- QV2M = (iy2 _ ]Vf2)/2M. 

Let us write the polarized structure functions as the sum of three contributions, 
namely 



gi{x,Q^ 



9f\x,Q 



'' + 9t''\x.Q') 
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(non-res.) 



x,Q') 



(15) 



(el.) 



where gf-'{x,Q'^), gr'-'{x,Q'^) and g. 



(res.) j 



(non—res.) 



X, Q"^) are the elastic, resonant and non- 



resonant contributions to gi{x,Q^), respectively, and i = 1,2. In Eq. (|TB|) possible in- 
terference terms between the resonant and non-resonant contributions are neglected, since 
they are well beyond the scope of our phenomenological fit. The Nachtmann moments (^) 
can then be written as the sum of three corresponding contributions, viz. 



M«(Q') = k«(Q 



(el.) 
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Mi'\Q' 
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3.1 Elastic term 

The contribution of the elastic process to the polarized structure functions can be expressed 
in terms of the electric and magnetic Sachs form factors as 

g\ '{x,Q } = d{x-l}-GM{Q ) — (17) 

Z i + T 

92 {x,Q ) = 6{x-1)-Gm{Q ) — (18) 

Z i + T 

where r = Q^/4Af^. Therefore, one gets 

L J 2 [ 1 + r 
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where ^ei = ^{x = 1) = 2/[l + Jl + 1/r]. For the explicit evaluation of Eq. (|T9D we will 
make use of the parameterization of Ref. [^ for the proton elastic form factors, assuming 
a 5% uncertainty. 

3.2 Non- resonant contribution 

The non-resonant terms gi"'~^'^^{x,Q'^) can be written in the following form 

g^r"'-\x,Q') = ^^-(x,Q^) + ^^/^(x,Q^) 

gt''—-\x,Q') = -g^-{x,Q')+g^^{x,Q') (20) 

where g^"{x^Q'^) is the contribution arising from the transverse asymmetry Ai(x, Q^) [pro- 
portional to (cri/2 — ^"3/2)], while g^'^{x,Q'^) is the LT contribution coming from the asym- 
metry A2{x,Q'^) [proportional to <Jlt], viz. 

, ,^, _ (non— res.) (non— res.) 

a^^ixO') = ^^ ^ "' '" 

y ^'^^ 47r2ae™ 1 + 4M2xVg2 2 

LT ( r)'i\ IVlri. i C^ (non— res.) (0^\ 

^ ^""'^ ^ ~ 47r2a,^l + 4M%VQ2 2Mx'^^^ ^ > 

Quite recently ||2^, a description of the transverse cross section difference (o"i/2 — cr3/2) 
in DIS kinematics for both proton and neutron has been obtained in terms of a parameter- 
ization inspired by the work of Ref. [Q. There, a simple Regge-type approach, based on 
one pomeron and one reggeon exchanges, has been shown to be phenomenologically success- 
ful in describing both the unpolarized photoproduction cross section and the unpolarized 



DIS data off the proton. The main result of Ref. [^ is that it is possible to parameterize 
smoothly the transition from the Regge behaviour, expected to be dominant at low values 
of Q^, to the partonic description valid at high values of Q^. Thus, one can try to use 
the same kind of parameterization to describe g^"{x^Q'^). That was the aim of Ref. [^ 



but unfortunately the explicit form of the parameterization adopted in ^^ fails just in the 
abovementioned smooth transition, because it does not possess a well defined Bjorken limit. 
Therefore, we have developed a new interpolation formula, properly inspired to the work of 
Ref. ^ 
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where W-,^ = M + m.,^ is the pion production threshold and t is a parameter aimed at 
describing the logarithmic scaling violations in the DIS regime, which we define following 



Ref. 30 as 



t = In 



. In (gg/A^) 



(23) 



In Eq. 



the parameter Q\ describes the transition from the expected dominance of the 



Regge behaviour at Q < Qj^ to the partonic regime at g >> Qj^. Indeed, on one hand, 
for g^ < g^ and at large W^ (i.e., low x) the r.h.s of Eq. (p2|) becomes proportional to 
J^jCijW^"'^ as expected from the Regge approach. On the other hand, when the Bjorken 
limit (fixed x and high g^) is considered, one gets g^°'{x, Q"^) oc Z^jLi djit) x~"j(*)(1 — x)'^^'^*^ 
Finally, the quantities a^, a^ 
on t, namely 



and j3j are parameters which are assumed to depend linearly 
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The parameters appearing in Eqs. (|22| - p3D have been determined by fitting existing 
measurements from Refs. |^T|-|^ on the asymmetry Ai(a;, g^) in the DIS kinematics, i.e. 
outside the resonance regions (W > 2 GeV), through the relation 



Ai{x,Q^)-*Dis 1 + 



AM^x^\ ^^"(x,g2 



g2 / Fi(x,g2) 



(25) 



where Fi(x, g^) is the unpolarized proton structure function calculated adopting the inter- 
polation formula developed in Ref. |]T|. Following Ref. |3^ 
Eq. (p^) fixing the QCD parameter A at the value A 



we have considered A^ = 2 in 
0.250 GeV. It turns out that our 
fitting procedure of DIS data is not very sensitive to the precise value of the parameter 
Wt, since it appears in Eq. (|22|) only in the combination Q^ + W^. Thus, the value of the 



9 



parameter Wt is relevant only at low Q^ and, indeed, in the next subsection it will be fixed 



by requiring the reproduction of the experimental value of the proton DHG sum rule [|15 
thus, we anticipate here its final value equal to Wt = 0.475 GeV. To sum up, we have used 
14 parameters against a total number of experimental points equal to 209, obtaining for the 
X^ variable divided by the number of d.o.f. the minimum value of 0.66. The explicit values 
of our parameters are reported in Table 1, where it can be seen that the transition point 
from the Regge behaviour to the partonic regime occurs around Q^ ~ Q|j ~ 4^5 {GeV/cY. 



Table 1. Values of the parameters appearing in Eqs. 
described in the text. 
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The quality of our fit is shown in Fig. 2. It can be seen that the differences between 
our fit and the data are approximately distributed as a gaussian-like distribution, and more- 
over they almost do not exceed the statistical + systematic errors of the data (added in 
quadrature). The uncertainty on our parameterization of g^'^ generated by the fitting proce- 
dure has been estimated through the uncertainties obtained for the values of the parameters 
reported in Table 1 from our least-x^ procedure. Using different bins in the variable x we 
have found that the uncertainty on our g^'^ can be approximated by the following simple 
formula: (1 — 0.5 ■ x/x-,,) ■ 15%, where x,r = Q"^ /{Q'^ + W^ — M^) is the pion threshold in x. 



We point out that in Eq. (|22D g is assumed to behave in the Bjorken limit as 
a power of x at low values of x. Since there is no strong argument in favor of such an 
assumption (cf., e.g., the discussion of the limit x ^ in Refs. |l^, 0), Eq. (|2^) has to 
be considered as a simple approximation valid in a limited x-range. In this respect we have 
already observed from Fig. 1 that existing data are scarce below x ~ 0.02; thus, we consider 
X > 0.02 as the x-range of applicability of our parameterization (P^, which is clearly well 
enough for our main purpose to study leading and higher twist effects at large-x. 

The contribution of g^'^{x,Q^) to Eq. (^Of) is marginal in DIS kinematics and, 
moreover, the whole effect of g2{x, Q^) in the Nachtmann moments (^ is power suppressed; 
therefore, we do not need a very refined interpolation formula for g2{x,Q'^)- In this respect 
the analyses of the SMC @, EU3 |21| and E155 |32| DIS data on A2(x, Q^) suggest that 
g2{x,Q^) is consistent with the Wandzura-Wilczek relation [^]. Thus, we impose a simple- 
minded generalization of the latter as a constraint on our parameterization for g^'^{x,Q'^); 
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Figure 2. (a) Distribution of the differences between our fit (^) and existing data |2l[| -[27| on the 
asymmetry Ai{x,Q'^) in DIS kinematics only (W > 1 GeV). (b) Ratio of the differences in (a) 
with the statistical + systematic errors of the data (added in quadrature) versus the final invariant 
mass W. 



we simply assume that 
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where Ft^r = v 1 — W^/W^ is a threshold factor. Thus, our interpolation model for the 
non-resonant contributions to gi{x,Q^) reads as 
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The corresponding predictions for A2{x,Q'^) = 2Mx ■ [gi{x,Q'^) + g2{x,Q'^)]/[Q ■ Fi{x,Q'^)] 
compare positively against available DIS data, as shown in Fig. 3. 



3.3 Resonant contribution and the DHG sum rule 



In the nucleon-resonance regions (W < 2 GeV) we follow the approach of Refs. [34, 35] by 
adopting a simple Breit-Wigner shape to describe the PF-dependence of the contribution of 
an isolated resonance R, while its Q^-dependence can be conveniently expressed in terms of 
transverse helicity amplitudes A^,2 and A^^m as well as the longitudinal helicity amplitude 
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Figure 3. Values of the asymmetry A2{x, (Q^)) for the proton versus the Bjorken variable x in 
case of DIS kinematics {W > 2 GeV). The open dots, squares and diamonds correspond to the 
data of the SMC [|2|, EIA2, |2|] and Elbh ^ experiments, respectively. The full dots are our 
predictions based on Eqs. (p7|- pq). The values of (Q^) are those corresponding to the kinematics 
of the various experiments for each value of x. 
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where Mr and T r = Tr{W = Mr) are the mass and the width of the resonance R, 



respectively, and eventually qw = \/{W^ + M^ - mD/AW^ - M^, qr = q^r(W = Mr) and 
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Qcm = yQ"^ + iW'^ ~ ^^ ~ Q'^Y / AW^ . The parameters (.r and Xr have the same meaning 
as in Ref. g. 



We point out that the VT-shape encoded in Eqs. ( PT| - |55D is inspired by the results 



of Ref. |24[, and therefore it differs from the one adopted in Ref. |3^, mainly because the 
width parameter Tr is assumed to be a constant independent on W in [^] at variance with 
Eq. (|33|). The use of the latter however allows to achieve a proper reproduction of the 
asymmetric shape around the resonance bumps due to PF-dependence of the available phase 
space for resonance decays. Such a dependence produces an important shift of the location 
of the resonance peaks with respect ioW = Mr, which is neglected in Ref. [^ . 

To develop our interpolation formula in the resonance regions we have included in 
Eqs. ( p^p^ ) all the "four-star" resonances of the PDG having a mass Mr < 2 GeV 



and a total transverse amplitude i/|^f/2P + 1^3/2 P larger than 0.050 GeV ^'^ at the photon 
point. Following Ref. [^ we have represented the Q^ behavior of the transverse helicity 
amplitudes in the following form 



|,K,,,,| . /i±4mc,e-«'. (34) 

while for the LT cross section we have introduced the parameter Af^ defined as 

^LT = -jR- (35) 

^1/2 

SO that crfy = A^rp ■ {Q/^/2qcm)o'^/2- Note that in Eq. ( p^ the parameter Cr represents 



just the total transverse amplitude y |^f/2p + 1^3/2^ ^^ ^^^ photon point. Finally, as for 
the background under the resonances we adopt directly Eq. (^) without any adjustment of 
the parameters, thanks to the fact that the latter smoothly behaves for W < 2 GeV. 

Firstly we have taken the values of the widths T^^ , of the amplitudes Cr and of 
the asymmetry A^ at the photon point from PDG [^. We have then calculated the 
asymmetry of the transverse cross section (T3/2 — o'i/2 as a function of the photon energy 
E^ = [W"^ — M^)/2M, and compared our results with the recent data from Mainz [Q. It 



turns out that an adjustment of the widths r^|^232) ^^^ -'^1)13(1520) ^^ ^^^^ ^^ °^ ^^^ amplitude 

Cdi3(i52o) is needed. Instead of their PDG values (put in parenthesis) we adopt: r^/^232) ~ 

0.10 (0.12) GeV, rgj^(^520) = 0-20 (0.13) GeV and Ci5,3(i52o) = 0.200 (0.170) GeV-^/'^. Fi- 
nally, as anticipated in the previous subsection, the value of the parameter Wt appearing in 
Eq. ( p2D is fixed at the value Wt = 0.475 GeV in order to reproduce the experimental value 
of proton DHG sum rule |15|, viz. 



Cr3/2- 0-1/2 27r2ag„K2 

dE^ — '-— — = — ~ 204.5 (xbarn (36) 

where E^^ is the pion threshold in terms of the photon energy E^ and k, is the (proton) 
anomalous magnetic moment. Our final results at Q^ = are reported in Fig. 4 and posi- 
tively compared with the Mainz data |T^. It can clearly be seen that our parameterization 
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of the non-resonant term ( pT] ) obtained from fitting DIS data, can be safely extended down 
to the photon point to describe a smooth background under the resonance bumps. 
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Figure 4. Asymmetry of the proton transverse cross sections, [(T3/2 — o"i/2], versus the photon 
energy E^. Fuh dots are the experimental data from Ref. |14|. The shaded area is our prediction, 
as explained in the text, while the dashed line is our non-resonant contribution. 

In Table 2 we have reported the contribution of various integration regions over the 
photon energy E'^ to the proton DHG sum rule ( p6D and to the forward spin polarizability 
7o, defined as 



7o 



47r2 



- ^3/2 



0"l/2 



(37) 



It can be seen that an important contribution to both the DHG sum rule and the forward 
spin polarizability comes from the experimentally unobserved region below E^ = 0.2 GeV 
in overall agreement with the prediction of the Unitarity Isobar Model of Ref. ||36|, while 
the contribution of our high-energy tail above E^ = l.Q GeV to the DHG sum rule is just 
half of the prediction of Ref. [^. Note that our parameterization of g^'^ (p^ ) implies that 
at Q^ = the transverse cross section asymmetry [(T3/2 — 0-1/2] 



behaves as E^ at high 



photon energies. The negative valus obtained for the parameters aj (see Table 1) largely 
ensure that the DHG integral (^) is convergent at high energies and does not require 
subtractions. 

The Q^ dependence of Af and the values of the parameter B^i and Af^j, have been 
estimated using available electroproduction data from Ref. [0. The uncertainty on our 
parameterization in the resonance regions has been estimated by assigning to the slope pa- 
rameter Bn an overall 25% uncertainty and to the parameter Gr the uncertainty AGr arising 
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Table 2. Contributions of various integration regions over the photon energy E^ (in GeV) to the 
proton DHG sum rule (B^) and to the forward spin polarizability 70 (pTJ), calculated using our 
interpolation formula. 



integration region 


DHG ifibarn) 


70 (10-6 fm*) 


E^ < 0.2 


-37 ±5 


133 ±17 


0.2 <E^< 0.8 


207 ±20 


-173 ±10 


0.8 <E^< 1.6 


47 ±9 


-5±1 


E^ > 1.6 


-13 ±2 


<0.1 


total 


204 ± 23 


-45 ± 20 



from the reported PDG uncertainties on the transverse helicity amplitudes at the photon 
point. Table 3 collects the values of all the relevant parameters necessary for evaluating Eqs. 
( pT| - |33|) . We have to mention that: i) the values adopted for the parameters Af;rp fully sat- 
isfy the bounds imposed by the existing data on the unpolarized longitudinal to transverse 
ratio (see Ref. |0), and ii) our parameterization of the full structure function g2{x, Q^) is 
consistent with the Burkhardt-Cottingham sum rule pOf, stating that /q dx g2{x,Q'^) = 0, 
within 2cr standard deviations at all values of Q^. 



Table 3. Values of the parameters appearing in Eqs. (31-33) and obtained as described in the 
text. Resonance masses are in MeV, Cr in GeV^^'^, AC/j in %, Br in GeV~^, ^r , Xr and Q 
in GeV. 



resonance 


Cr 


^Cr 


Af 


Br 


r(o) 


aR 


iR 


Xr 


A(1232) 


0.290 


5 


-0.56 


0.7 


0.10 


-0.1 


1 


0.16 


Pii(1440) 


0.065 


5 


1.0 


1.6 


0.30 


-0.2 


1 


0.35 


L)i3(1520) 


0.200 


10 


-^ _ g0.67-0.65g 


0.8 


0.20 


0.2 


2 


0.35 


5ii(1535) 


0.090 


35 


1.0 


0.6 


0.15 


0.2 


1 


0.35 


5i*i(1650) 


0.055 


30 


1.0 


1.0 


0.15 


-0.3 


1 


0.35 


Di5±Fi5(1680) 


0.150 


15 


-^ _ „0.69-0.98g 


0.6 


0.13 


-0.3 


3 


0.35 


/^33(1700) 


0.135 


30 


0.2 


1.0 


0.30 


-0.3 


2 


0.35 


^35(1905) 


0.055 


60 


-0.50 


0.6 


0.25 


-0.3 


3 


0.35 


^37(1950) 


0.125 


20 


-0.24 


0.6 


0.25 


-0.3 


3 


0.35 



The quality of our interpolation formula in the resonance electroproduction regions 
is illustrated in Fig. 5, where our predictions are positively compared with the two sets of 
existing data from EU3 at Q^ ^ 0.5 and ~ 1.2 (GeV/c)^. 



Before going on with the calculation of the Nachtmann moments (^, our prediction 
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Figure 5. Comparison of the proton transverse asymmetry A^(W,Q'^) obtained through our 
interpolation formula (full dots) with the existing data from Ref. |2^] (open dots) in the resonance 
regions. The average value of Q^ is ~ 0.5 and ~ 1.2 [GeV/c)'^ in (a) and (b), respectively. 



for the generalized DHG integral, defined as 

2M2 



I DHG 



iQ' 



Q' 



dx gi{x, Q^ 



(38) 



is presented in Fig. 6 and positively compared with available experimental data. Note that 
in Eq. (pSf ) the upper limit of integration excludes the elastic contribution and therefore 
the generalized DHG integral Idhg{Q^) cannot be analyzed in terms of twists, because the 
OPE is fully inclusive. 

As it is well known, at the photon point the (proton) generalized DHG integral 



iDHciQ^ = 0) 



M2 



STT^ttem "'^^ 



°°ciE"^/^~"^/^ 



E^ 



K 

T 



(39) 



is negative (~ —0.80), while it becomes positive in the DIS regime. Therefore it should 
cross zero at some value of Q"^ and according to the predictions of our parameterization (see 
Fig. 6) the zero-crossing point is expected to occur at Q"^ = Q\ = 0.16 ± 0.04 {GeV/cy. 



Such a value is significantly below the prediction Qx ~ 0.8 (GeV/c)'^ found in Ref. pO| , 
where the dominance of the N — A(1232) transition in Eq. (|38D at low Q^ was assumed 
and a vector meson dominance picture was adopted for the non-resonant background. The 
difference is due to the significative contribution of non-resonant processes present at low Q^ 
in our parameterization (see dashed curve in Fig. 6) at variance with the assumptions made 
in Ref. [|^. Our zero-crossing point is below the finding Qx ~ 0.3 (GeV/c)'^ obtained in Ref. 
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pT| within the constituent quark model and characterized by the so-called dominance of low- 
lying resonances at low Q'^ . Our result is quite close to the prediction Q\ ~ 0.2 (GeV/c)^ 
of Ref. [|2|, where the origin of the zero-crossing point is traced back to the strong Q^- 
dependence of the first moment of the inelastic part of the structure function g2{x,Q'^) 
implied by the Burkhardt-Cottingham sum rule [p9|. Therefore, our finding, if confirmed by 



a direct measurement of the Q^-behavior of the generalized DHG integral, should provide 
an important constraint on hadronic models and on the physics of nucleon resonances. 
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Figure 6. Generalized DHG integral (|3q ) for the proton versus Q^ 

dots and triangles are the experimental results of Refs. |^, 22, |2^, 

prediction based on our parameterization of g\{x,Q'^), while the dashed line corresponds to the 

non-resonant contribution only. The arrows indicates the location of the value of DHG sum rule 

(~ —0.80) at the photon point. 

The Nachtmann moments M^^Q"^) [see Eq. (^)] can now be evaluated through our 
interpolation formulae for the polarized structure functions gi{x,Q'^) and g2{x,Q'^), which 
we stress have the following range of applicability: x > 0.02 and Q^ < 50 [GeV/cY. The 
separate contributions of the elastic peak (|r^-|TB|), the resonances ( PU| - P^ and the inelastic 
contribution [i.e. the sum of the resonant and non-resonant parts (p7|-p^)] are shown in 
Fig. 7 for 0.1 < Q^ {GeV/c)"^ < 50. It can be seen that the elastic term is dominant for 
Q^ ^ ^/2 {GeV/cy, while the inelastic contribution exhausts the Nachtmann moment for 

H), the 



Q >, n (GeV/c) . Note that, at variance with the generalized DHG integral 



first moment MJ {Q^) 

contribution (|19|), as firstly pointed out in Ref. 



does not become negative at low values of Q^ because of the elastic 

|. Moreover, the inelastic parts of the 
Nachtmann moments shown in Fig 7(b-d) change their sign at Q^ ~ 0.45 ± 0.10, 
0.15, 1.1 ± 0.2 {GeV/cf for n = 3, 5, 7, respectively. 
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Figure 7. Proton Nachtmann moments MA (Q^) [see Eq. (|3|)] versus Q^ for n = 1 (a), n = 3 
(b), n = 5 (c) and n = 7 (d). Open diamonds and triangles correspond to the contributions 



of the elastic peak (IT-p^) and the resonances ( p9[]3^ ), respectively. Open dots are the inelastic 
contribution, i.e. the sum of the resonant and non-resonant parts (27-|2^). Open squares represent 
the full Nachtmann moment given by the sum of the elastic and inelastic parts. 



Finally, we want to mention that quite recently [^ preliminary photoabsorption 
data above E^ ~ 0.8 GeV have become available from a GDH experiment at ELS A. Such 
data cover a large kinematical region where the Fi5(1680) nucleon resonance contributes 
predominantly. It appears that our fit, based on the PDG values for the helicity amplitudes 
at the photon point, overestimates the preliminary data. The agreement can be recovered 
simply by imposing a ~ 25% reduction of the strength parameter (^^^^(leso) and a ~ 10% 
reduction of the parameter Wt in order to keep the GDH sum rule fulfilled. We do not 
include these modifications in our present parameterization, leaving this issue to be fixed 
when final data from ELS A will be available. Nevertheless, we have checked that the above- 
mentioned changes do not modify within the quoted uncertainties both the zero-crossing 
point of the generalized GDH integral and the twist analysis of the next two Sections. 



3.4 Bloom-Gilman local duality 



An important feature of the results shown in Fig. 7 is that the resonant contribution is 
negative for Q^ ~ few {GeV jc)^ . This is mainly due to the well established fact [^ that the 
proton transverse asymmetry J^^ in the regions of the A(1232)-resonance electroproduction is 
negative up to Q^ ~ 3-i-4 {GeV jc)^ . Thus, for Q^ ~ few [GeV jc)^ the resonant contribution 
to the polarized proton structure function is opposite in sign with respect to the unpolarized 
case (see Ref. 0]). In this respect we want ot remind that the concept of parton-hadron local 
duality makes an important (and not yet fully understood) connection between the physics 
in the nucleon-resonance and DIS regions. Indeed, the parton-hadron local duality, observed 
empirically |T^ by Bloom and Gilman in the unpolarized transverse structure function of 
the proton, states that the smooth scaling curve measured in the DIS region at high Q^ 
represents an average over the resonance bumps seen in the same x region at low Q^ . More 
precisely, it occurs a precocious scaling of the average of the F^i^, Q^) data in the resonance 
regions to the DIS structure function -Ff(.^), at corresponding values of the improved scaling 
variable ^ p5[. 
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Figure 8. Polarized proton structure function g^H^Q"^) versus the Nachtmann variable ^ for 
various values of Q^. In (a)-(d) the shaded areas are pseudo-data generated through our interpo- 
lation formula at Q^ = 0.5,1,2 and 5 {GeV/c)^, respectively. The solid line is the result of our 
parameterization evaluated at Q^ = 20 {GeV/c)^ . 
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It is therefore legitimate to ask ourselves whether the BG local duality holds as well in 
the polarized case. To this end we have generated pseudo-data in the resonance regions and 
in the DIS regime via our interpolation formula for gi{x,Q'^). Our results are reported in 
Fig. 8, where it can clearly be seen that: i) at values of Q^ as low as ~ 0.5 {GeV/cY there is 
no evidence at all of an occurrence of the BG local duality, as in the case of the unpolarized 
transverse structure function of the proton (see Ref. 0)), and ii) in the kinematical regions 
where the A(1232) resonance is prominently produced, the BG local duality breaks down at 
least for Q^ up to few {GeV/cY, while in the higher resonance regions for Q^ > 1 [GeV/c)'^ it 
is not excluded by our parameterization. Note that in the unpolarized case the onset of the 
BG local duality occurs at Q^ ~ l-=-2 (GeV/c)'^ HO' El' including also the A(1232) resonance 
regions [Q . It should be mentioned that the usefulness of the concept of local duality relies 
mainly on the possibility to address the DIS curve at large ^ through measurements at low 
Q^ in the resonance regions. It is therefore clear that the breakdown of the local duality 
in the region of the A(1232) resonance forbid us to get information from duality on the 
behavior of the scaling curve at the highest values of ^. 

4 NLO analysis of the polarized Nachtmann moments 

In this Section we present our power correction analysis of the polarized Nachtmann mo- 
ments @, adopting for the leading twist the NLO approximation. Following Refs. [|l], ^ a 
phenomenological ansatz is introduced for describing power corrections, viz. 



Mi'\Q') = Sf,l'\Q') + Sa^:^ 









6a^'^ 






X (6) 



(40) 



where the leading twist term 5ix^\Q'^) is given by Eqs. (|9|-|TTD, while the logarithmic pQCD 
evolution of the twist-4 (twist-6) contribution is accounted for by the term [as{Q^) / as{fi'^)Y'^^ 
([as(Q^)/a;s(/i^)]^'^i ) with an effective anomalous dimension 57} (57} ) and the parameter 
5ai {5a[ ) represents the overall strength of the twist-4 (twist-6) term at the renormaliza- 
tion scale /z^. For the latter we consider hereafter the value fx = 1 GeV/c and, for fixing 
the running of the coupling constant ^^(Q^), the updated PDG value as(M|) = 0.118 ||^ is 
adopted throughout this work. 

In Eq. (|40|) only twist-4 and twist-6 terms are included. In this respect we want 
to point out that the number of higher-twist terms to be considered is mainly governed by 
the Q^-range of the analysis. Indeed, as the latter is extended down to lower and lower 
values of Q^, more and more higher- twist terms are expected to contribute equally well. We 
anticipate here that for Q^ > 1 {GeV/c)"^: i) the inclusion of a twist-4 and a twist-6 term 
appears to work pretty well, as already found in case of the unpolarized moments |l|, 0, 
and ii) our least-x^ fitting procedure turns out to be not sufficiently sensitive for a precise 
determination of power corrections of order higher than the twist-6. Thus, just because of 
phenomenological findings we limit ourselves to consider only twist-4 and twist-6 terms in 
our analyses for Q^ > 1 {GeV/cY. 

20 



Eqs. 



Let us start with the twist analysis of the first Nachtmann moment MJ (<5^); from 
0|) and ([10|) one has 



r(i) 



Ml'\Q'' 



< e^ > 



2 



Ag^^^+ao(Q^ 



1 - 






X (4) 
2\1'57i ,,2 



/^ , X (6) 

0^ + ^"^ 






2M'57i 



as(/i^ 



(6) 




(41) 



The non-singlet moment Ag^"^ is taken fixed at the value Ag^"^ = 1.095, deduced from 
the experimental values of the triplet and octet axial coupling constants [see Eq. (]T2|)], 
with the latter obtained under the assumption of 5't/(3)-flavor symmetry. The unknown 
parameters in Eq. (^Tj) are the singlet axial charge ao(/i^) at the renormalization scale (or 
at any given value of Q'^) and the four higher-twist quantities 5a\ , ^7! , 5ai and 5'^i . 
Their values, reported in Table 4, have been determined by fitting the pseudo-data of Fig. 
7(a), adopting the least-x^ procedure in the Q^-range between 0.5 and 50 {GeV/cY. The 
twist decomposition of the Q^ behavior of the first moment M^ (Q^) is illustrated in Fig. 9. 



Table 4. Values of the parameters appearing in Eq. (|4l| ) obtained by a least-x^ procedure in the 
Q^-range from 0.5 to 50 {GeV/c)'^. The non-singlet moment is fixed at the value Aq = 1.095 
(see text). The last row reports the minimal value obtained for the x^ variable divided by the 
number of degrees of freedom. The errors on the parameters represent the uncertainty of the fitting 
procedure corresponding to one-unit increment of the X^ /^d.o.f. variable. 



ao(10 GeV^) 


5a^^^ 


X (4) 


5af) 


Sir' 


xVN,.,j. 


0.14 ±0.09 


0.038 ±0.012 


2.2 ±0.4 


-0.017 ±0.006 


1.9 ±0.6 


0.053 



From Table 4 and Fig. 9 it turns out that: 

• the total contribution of the higher twists is tiny for Q^ > 1 {GeV/cy, but it is 
comparable with the leading twist already at Q^ ~ 0.5 {GeV/cY. Since the first 
moment basically corresponds to the area under the structure function Qi (as it is the 
case of the second moment of the unpolarized structure function Ff), the dominance of 
the leading twist in Mf (Q^), occurring for Q^ >l {GeV/cf', reflects only the concept 
of global duality and not that of local duality (cf. Ref. [H); 

• in our analysis, where the leading and the higher twists are simultaneously extracted, 
the singlet axial charge (in the AB scheme) is determined to be ao(10 GeV'^) = 0.14 ± 
0.09 (see Table 4), which nicely agrees with many recent estimates appeared in the 
literature, like e.g. ao(10 GeV^) = O.lO+g^I from Ref. |T6[ and ao(10 GeV^) = 
0.24 ± 0.07(stat) ± 0.19(syst) from Ref. lUKb). As a consistency check, we have 
limited our analysis to the high Q^-range from 10 to 50 {GeV/c)"^ including only the 
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Figure 9. Twist analysis of the first proton Nachtmann moment M} (Q^). The solid line is the 
result of Eq. (41) fitted by the least-x^ procedure to our pseudo-data of Fig. 7(a) (open squares) 
in the Q^-range between 0.5 and 50 {GeV/c)^. The dotted and dashed lines correspond to the 
contributions of the leading and total higher twists, given by the sum of the twist-4 and twist-6 
terms in Eq. (|4l]), respectively. 



twist-2 contribution in Eq. (^), obtaining ao(10 GeV"^) = 0.18 ± 0.09. Thus, possible 
higher-twist effects on the extraction of ao are small and well within the uncertainties 
of the fitting procedure. We quote ao(10 GeV^) = 0.16±0.09 as our final determination 
of the singlet axial charge in the AB scheme. Our value of oq is therefore significantly 
below the naive quark- model expectation (i.e. compatible with the well known "proton 
spin crisis"), but it does not exclude completely a singlet axial charge as large as ~ 0.25. 

As explained in Section 2 [see Eq. (P)], for higher-order moments {n > 3) we make 
use of the following twist expansion 



r(l) 



Ml^UQ' 



SAi-^ 




1 + ^5C(^) 

27r 






+ 6a^S 



as(/i2 



Sji'' .2 



Q 







(42) 



The five parameters SA'^jj^), 5a^\ ^In^ ^'^n^ ^^^ ^In^ ^^^ simultaneously determined by 
the least-x^ procedure apphed to the pseudo-data of Fig. 7 in the Q^-range from 1 to 
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50 {GeV /cY; the obtained values are reported in Table 5, while the twist decomposition of 
mS3(Q2) is illustrated in Fig. 10. 



Table 5. Values of the parameters appearing in Eq. ( p^ at /x = 1 {GeV/c), obtained by a least-x^ 
procedure in the Q^-range from 1 to 50 {GeV/c)'^. The errors on the parameters represent the 
uncertainty of the fitting procedure corresponding to one-unit increment of the X^/^d.o.f. variable. 



parm. 


n = 3 


n = 5 


n = l 


n = 9 


M„ 


0.133 ±0.019 


0.0297 ±0.0044 


0.0098 ±0.0013 


0.00343 ±0.00037 


5aif) 


0.015 ±0.005 


0.018 ±0.005 


0.026 ±0.003 


0.032 ±0.004 


Slk'^ 


1.7±0.6 


2.4 ±0.9 


3.7 ±0.5 


4.6 ±0.5 


^«f 


-0.006 ±0.002 


-0.014 ±0.005 


-0.024 ±0.004 


-0.031 ±0.005 


Sll^^ 


1.7 ±0.6 


1.7 ±0.5 


2.6 ±0.5 


3.3 ±0.5 


x'/Nd.o.f. 


0.14 


0.39 


0.60 


0.82 



Our main results can be summarized as follows: 

the twist-2 term, extracted from our proton pseudo-data simultaneously with the twist- 
4 and twist-6 contributions, differs only slightly from the predictions obtained using 
the polarized PDF set of Ref. |^ evolved at NLO (compare crosses and dotted 
lines in Fig. 10). A similar situation holds as well in case of the PDF set of Ref. 
0. Moreover, we have checked that, by limiting our analysis to the Q^-range from 
10 to 50 {GeV/cY and without including any higher-twist term, our extracted twist-2 
changes only within the fitting uncertainties reported in Table 5; 

the twist-4 and twist-6 contributions turn out to have opposite signs, making the total 
higher-twist contribution smaller than its individual terms, as we have already found 
in case of the unpolarized proton and deuteron structure functions (see Ref. |1|. As n 
increases (i.e. as x increases), the total higher- twist contribution increases and becomes 
comparable or even larger than the leading twist term for Q^ ~ few {GeV/cf' (compare 
dashed and dotted lines in Fig. 10); 

for n = 3 -h 9 the values of the effective anomalous dimensions 57^"^^ (and to a less 
extent also the values of 6'y^^) result to be larger than the values of the corresponding 
twist-2 anomalous dimensions (7^"^ ~0.6^1.2 forr2 = 3-r-9 



the uncertainties on the different twist contributions due to our least-x^ procedure are 
within ~ 40% for the twist-4 and twist-6 terms, while they are about 10 -h 15% for the 
leading twist (see Table 5); 

the twist expansion ( ^21) appears to work quite well for values of Q^ down to ~ 
1 {GeV/cf. 
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Figure 10. Twist analysis of the proton Nachtmann moments M„ {Q ) with n = 3 (a), n = 5 



(b), n = 7 (c) and n = 9 (d). The sohd line is the result of Eq. ( |4^ ) fitted by the least-x^ 
procedure to our pseudo-data (open squares). The dotted, dot-dashed and triple-dot-dashed lines 
are the separate contribution in Eq. (|4^ ) of the twist-2, twist-4 and twist-6 terms, respectively. The 
dashed line is the total higher-twist contribution, given by sum of the twist-4 and twist-6 terms. 



The crosses are the moments calculated using the NLO polarized PDF set of Ref. [47|, labeled 
valence scenario. 



5 Soft gluon resummation and higher twists 



So far the power corrections appearing in Eqs. ( |4l| - ^2D have been extracted from our pseudo- 
data assuming the NLO approximation for the leading twist, and therefore they represent 
the higher twists at NLO. Since our main aim is to get information on the dynamical 
power corrections generated by multiparton correlations, it is necessary to estimate the 
possible effects of higher orders of the perturbative series, which defines the twist-2 coefficient 
functions En2[fJ',o:s{Q'^)] appearing in Eq. (|]). Such a job has been already carried out in 
case of the unpolarized proton structure function F^x, Q^) in Refs. |]I[ and 0. In the former 
the high-order perturbative terms have been estimated through the ambiguities introduced 
by infrared renormalons, while in the latter the effects of soft gluon emission have been taken 
into account at large x via well-established soft-gluon resummation techniques 0. In both 
cases the comparison of the results obtained adopting the NLO approximation and those 



24 



including high-order corrections has clearly shown that the extraction of higher twists at 
large x is remarkably sensitive to soft gluon effects. Therefore, in this Section we address 
the issue of high-order effects by extending the calculations of Ref. [^ to the Nachtmann 
moments of the polarized proton structure function gi{x, Q^). 

Let us start from Eq. (^) [see also Eq. (|^)] where 6C!;^^ is the NLO part of the 
quark coefficient function. In the MS (or AB) scheme it reads explicitly as 



SCi^^ = CF\s,in) 






In n + 1 n'^ 



(43) 



where Cp = {N^ - l)/(2iV,) -^ 4/3 and Sk{n) = E"=i 1/j^ Note that SC^^^ differs from 
the corresponding quark coefficient function of the unpolarized case (cf., e.g., Eq. (4) of Ref. 
1^) only in the terms proportional to 1/n and l/in + 1) in the r.h.s. of Eq. (|45|). Thus, the 
SGR effects can be included in the polarized NS moments in a way completely similar to 
the unpolarized NS case already carried out in Ref. ||^. For completeness here below we 
report the explicit expressions of SGR effects in the polarized NS case. 

For large n (corresponding to the large-x region) the coefficient 5G^^ is logarithmically 
divergent. Indeed, since S'i(n) = ■Je + log(n) + 0(l/n) and S2{n) = 7r^/6 + 0(l/n) (with je 
being the Euler-Mascheroni constant), one gets 



SGi^^ = 6G]^>s + SCl^LOG + 0{l/n) 



(<?) 



where 



SG 



DIS 



TX 



5G, 



n,LOG 



Cf rE + ^lE~--- 



GfHti) ln(n) + 27£; + - 



(44) 

(45) 
(46) 



The physical origin of the logarithm and double logarithm terms in Eq. ( ^61) is the 
mismatch among the singularities generated in the quark coefficient function by the virtual 
gluon loops and the real gluon emissions, the latter being suppressed as the elastic peak 
corresponding to the threshold x = 1 is approached. In other words at large x the relevant 
scale is no more Q^, but Q^(l — x) |^0|| and the usual Altarelli-Parisi evolution equation 
should be accordingly modified [^. The presence of the abovementioned divergent terms 
at large n would spoil the perturbative nature of the NLO approximation (as well as of any 
fixed-order calculation) and therefore the effects of soft gluon emissions should be considered 
at all orders in the strong coupling constant a^. To this end one can take advantage of 
resummation techniques, which show that in moment space soft gluon effects exponentiate 
0, ^, Q. Thus, the moments of the leading twist including SGR effects acquire the 
following form 



r (1,SGR) 



{Q') = SA^ifi' 
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27r 

asm 



-C 



(g) 



DIS 

as{iJ?) 



gCnCQ^) ^ 



as{Q^ 

27r 



' r^{q) r^{>i) r^ii) 



in 



% 



l,NS 



A ATS 



(47) 



where the function Gn{Q^) is the key quantity of the soft gluon resummation and reads as 

(cf. mM) 



GniQ' 



1 ^"-1 

dz 

1- 



1 (1 



BMQ'il-z))] 



g2 



A[asiq')] 



(4J 



where A[a,] = Cpa,/7r + (:7pir(a,/7r)V2, 5[a,] = -3(:7^a,/27r with K = C^ (67/18 -7rV6)- 
10TrA^//9, Ca = Ne-^3 and Tr = 1/2. Exphcitly one has 



G„,(g2) = ln(n)Gi(A 

where A„ = /3oas(Q^)ln(n)/47r and 

4 



G2(A„) + 0Kln'=-^(n)] 



(49) 



G'i(A) 
G'2(A) 



Cf^tt[^+ (1-^)1^(1-^)] 



-C 



/3oA 

47E + 3 



/5o 



ln(l - A) - Cf^ [A + ln(l - A)] + 

Po 



C, 



4A 



A + ln(l-A) + -ln2(l-A) 



(50) 



.(<?) 



It is straightforward to check that in the hniit A„ << 1 one has GniQ"^) — ^ as{Q'^)C)^j^QQ/2T[, 
so that at NLO Eq. (^) reduces to Eq. (|^). Note that the function G2(A) is divergent 
for A — >■ 1; this means that at large n (i.e. large x) the soft gluon resummation cannot be 
extended to arbitrarily low values of Q^. Therefore, for a safe use of present SGR techniques 
we will work far from the above-mentioned divergencies by limiting our analyses of low-order 
moments {n < 9) to Q^ > i [GeV/cf Q. 

The results obtained for the ratio of the quark coefficient function calculated within 
the SGR technique and at NLO, namely 



TniQ' 



1 + a.(g^)Ca,/2vr) e^-(Q^) + a^m (cj^) - C^^]s - ^^0^) /27r 



l + «,(g2)CA''727r 



(51) 



are reported in Fig. 11 at ag^M'^) = 0.118. It can be seen that soft gluon effects are quite 
small for the ffist and the third moments, while for n > 5 they increase significantly as n 
increases, particularly for Q^ ~ few [GeV/cY. Note that for n = 1 one has ri(Q^) = 1 at 
any values of Q^ and as{M'^)\ therefore, no effect from soft gluons can occur in the ffist 
moment m[^\Q'^). 

''We want to stress that, as in Refs. |0, 0], our main aim is not to perform a full calculation of perturbative 
corrections, but to check whether the NLO approximation can provide a safe extraction of higher twists. 
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Figure 11. Values of the ratio rn{Q'^) [see Eq. (|5l|)] of the quark coefficient function within the 
SGR technique and at NLO. The dotted, dot-dashed, triple-dot-dashed, dashed and solid lines 
correspond to n = 1,3,5,7 and 9, respectively. Note that ri{Q'^) = 1 and r3((5^) < 1. 



Thus, for n > 3 our pseudo-data of Fig. 7 have been analyzed in the Q^-range from 
1 to 50 {GeV/cY including soft-gluon effects, viz 



M«(Q^) = 5^^"^^HQ^) + 5«lf) 






<57i^> 









-Si's' 



k) <-' 



where the leading twist term 5jji\l'^'^^\Q'^) is given by Eq. (0). The five parameters 
5An{iji^)-, 5a^n''^ ^ln\ ^^n^ ^^^ ^7^^^ are simultaneously determined by the least-x^ procedure; 
the obtained values are reported in Table 6, while the twist decomposition of M„>3((5^) is 
illustrated in Fig. 12. It can be seen that the twist-2 parameters 5 An (see Table 6) almost 
coincide with the corresponding NLO quantities 5 An (see Table 5) within the uncertainties 
generated by our least-x^ procedure. As in the case of the NLO analysis, the twist-4 and 
twist-6 contributions are again well defined and with opposite signs, making the total higher- 
twist contribution smaller than its individual terms. The strengths 5a!^'> and daf^n^ (see Table 
6) differ remarkably from the corresponding NLO quantities Sa^n'' ^^^ 5a^^-* (see Table 5), 
while the effective anomalous dimensions 5'^^n^ and ^7^^^ do not change significantly. We 
point out that at large n (i.e. at large x), the total higher-twist contribution is reduced by 
soft gluon effects, but it is still a significant fraction of the leading twist term for Q^ ~ few 
{GeV /cY (compare dashed and dotted lines in Fig. 12). 

The comparison of our twist analyses at NLO and within the SGR technique is 
reported in Fig. 13, where it can be seen that, except for the third moment (which is the 
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Table 6. Values of the parameters appearing in Eq. (52) at /x = 1 {GeV/c), obtained by a least-x^ 
procedure in the Q^-range from 1 to 50 {GeV/c)^. The errors are as in Table 5. 



parm. 


n = 3 


n = 5 


n = 7 


n = 9 


5An 


0.135 ±0.025 


0.0297 ±0.0050 


0.0092 ±0.0016 


0.00358 ±0.00064 


5a^^^ 


0.015 ±0.005 


0.0080 ±0.0025 


0.0067 ±0.0008 


0.0080 ±0.0009 


^7i^^ 


1.9 ±0.7 


1.9 ±0.6 


2.4 ±0.9 


4.0 ±1.5 


H^) 


-0.0062 ±0.0022 


-0.0053 ±0.0019 


-0.0073 ±0.0007 


-0.0011 ±0.0002 


Hf^ 


1.9 ±0.8 


3.0 ±1.1 


3.5 ±1.2 


3.7±1.3 


xVN,.o.f. 


0.13 


0.29 


0.46 


0.67 



only one characterized by r^{Q'^) < 1), the contribution of the twist-2 is enhanced by soft 
gluon effects, while the total higher-twist term decreases significantly after the resummation 
of soft gluons. Thus, as already observed in the unpolarized case in Ref. @, also in the 
polarized one it is mandatory to go beyond the NLO approximation and to include soft 
gluon effects in order to achieve a safer extraction of higher twists at large x, particularly 
for Q2 _ few [GeV/cf. 

In Fig. 14 we have compared the twist decomposition of the polarized Nachtmann 
moments M!f^\Q'^) (with n = 3,5,7) obtained in this work and the corresponding decom- 
position of the unpolarized (transverse) Nachtmann moments M!^\Q'^) (with n = 4,6,8) 
obtained in Ref. adopting the same SGR technique. It can clearly be seen that our ex- 
tracted higher-twist contribution appears to be a larger fraction of the leading twist in case of 
the polarized moments. This findings suggests that spin- dependent multiparton correlations 
may have more impact than spin-independent ones. 

6 Conclusions 



In conclusion, we have extend the twist analysis made in Refs. [|l| and [^] to the case of the 
polarized proton structure function gi{x, Q^). Within the framework of the Operator Product 
Expansion we have used the Nachtmann moments in order to disentangle the kinematical 
target-mass corrections from the dynamical higher-twist effects related to correlations among 
partons. Since the evaluation of the Nachtmann moments requires the knowledge of the 
polarized structure functions in the whole x-range for fixed values of Q^, we have developed 
a new parameterization of gKxjQ"^), which describes the DIS proton data up to Q^ ~ 
50 {GeV/cy and includes a phenomenological Breit-Wigner ansatz able to reproduce the 
existing electroproduction data in the proton-resonance regions. Our interpolation formula 
for gi{x, Q"^) has been successfully extended down to the photon point, showing that it 



can nicely reproduce the very recent data |IJ] on the energy dependence of the asymmetry 
of the transverse photoproduction cross section as well as the experimental value of the 
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Figure 12. Twist analysis of the proton Nachtmann moments M„ {Q'^) with n = 3 (a), n = 5 



(b), n = 7 (c) and n = 9 (d). The sohd line is the result of Eq. ( [5^ ) fitted by the least-x^ 
procedure to our pseudo-data (open squares). The dotted, dot-dashed and triple-dot-dashed lines 
are the separate contribution in Eq. (p^ ) of the twist-2, twist-4 and twist-6 terms, respectively. The 
dashed line is the total higher-twist contribution, given by sum of the twist-4 and twist-6 terms. 



proton Drell-Hearn-Gerasimov sum rule. According to our parameterization of gi{x, Q"^) the 
generalized Drell-Hearn-Gerasimov sum rule is predicted to have a zero-crossing point at 
Q2 = 0.16 ±0.04 {GeV/cf. 

Low-order polarized Nachtmann moments have been evaluated and their Q^ behavior 
has been investigated in terms of leading and higher twists for Q^ > 1 {GeV/c)"^. In our 
analyses the leading twist is extracted simultaneously with a phenomenological higher-twist 
term from our pseudo-data. We have reported results obtained both at NLO and adopting 
the same soft-gluon resummation technique applied to the analyses of the unpolarized data 
made in Ref. in order to take into account the effects of higher orders of the perturbative 
series. As far as the first moment is concerned, the effects of higher twists are found to be 
quite small for Q^ > 1 (GeV/c)'^. Moreover, the singlet axial charge is determined to be 
ao[10 [GeV/cy] = 0.16 ± 0.09; our extracted value is significantly below the naive quark- 
model expectation (i.e., compatible with the well-known "proton spin crisis"), but it does 
not exclude completely a value of the singlet axial charge as large as ~ 0.25. In case of 
higher order moments, which are more sensitive to the large- a; region, higher-twist effects 
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are significantly reduced by the introduction of soft gluon contributions, but they are still 
relevant for Q^ ~ few {GeV/cY at variance with the case of the unpolarized transverse 
structure function of the proton examined in Ref. Q]. Our finding suggests that spin- 
dependent correlations among partons may have more impact than spin-independent ones. 
As a byproduct, it has also been shown that the Bloom-Gilman local duality is strongly 
violated in the region of polarized electroproduction of the A (1232) resonance. 

The analysis of deuteron data, aimed at the determination of the flavor dependence 
of leading and higher twists, is in progress and the results will be presented elsewhere. 

The results presented in this work and those already reported in Refs. [Q and P 
demonstrate that power correction terms can be safely extracted from a proper analysis of 
inclusive data. In particular, our extracted values of the effective strengths of the twist-4 and 
twist-6 terms may be directly compared with theoretical results obtained from first-principle 
calculations (like lattice QCD simulations) or from models of the nucleon structure. In this 
way the non-perturbative regime of QCD may be tested. 

Let us finally stress that our present analysis is mainly limited by the use of phe- 
nomenological fits of existing data (i.e., pseudo-data)^ which are required for interpolating 
smoothly the nucleon structure functions in the whole x-range for fixed values of Q^. There- 
fore, polarized inclusive data with better quality at x > 0.5 and Q^ < 10 {GeV/cY, which 
may be collected at planned facilities like, e.g., J LAB @ 12 GeV, could greatly help to 
improve our understanding of the non-perturbative structure of the nucleon. Finally, we 
want to point out that, since in inclusive polarized and unpolarized processes multiparton 
correlations appear to generate power-like terms with opposite signs, semi-inclusive or ex- 
clusive experiments might offer the possibility to achieve a better sensitivity to individual 
non-perturbative power corrections. 
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Figure 13. Comparison of the twist analyses of the proton polarized Nachtmann moments 
^iS-siQ^) at NLO [see Eq. (||)] and within the SGR technique [see Eq. (||)] for n = 3 (a,b), 
n = 5 (c,d) and n = 7 (e,f). Dotted lines are the twist-2 term, while dashed lines correspond to the 
total higher-twist contribution, given by sum of the twist-4 and twist-6 terms. The shaded areas 
help to identify the separate contributions of the leading and higher twists. 



34 



o 



o 



o 



0.03 



0.02 



0.01 



0.00 



— I 1 1 — I — I — r-r-| 1 — 

(a) pol.: n = 3 







10 

Q^ {GeV/cf 



0.010 



0.005 i/. 



0.000 



0.004 



0.002 



0.000 



1 1 1 — r 



-I 1 1 — I — I — I I I I — 

(c) pol.: n = 5 




Q^ {Gey/cf 



T 1 1 — r 



— I 1 1 — I — I — I I I I — 

(e) pol.: n = 7 




J I I L 



Q^ {GeV/cf 



O 



O 



0.04 



0.02 



0.00 



0.02 



(b) unpol. 







10 



Q^ {Gey/cf 



O 0.01 



0.00 



(d) unpoL: n = 6 




-\ \ _ 
y y y .* , 

\ \ X \ ■ 
y y y y y y" -. . 

\ \ ■-. \ \ ■-. \ ' ^ • ^ , 

\ \ \ \ \ \ \ \ \ \ \ \ --^ 
y y y^-v ^^ ^ y ^ ^ y ^ ^ y ^ ^ y ^ -r .- .- 






Q^ {Gey/cf 



10 



0.008 



0.004 



0.000 




\ X 

1^ y y 

\ X \ \- 

y y y y y" , 

\ \ \ \ \"p,i 
yyyyyyy 

\x\\x\\x" 
yyyyyyyyyyy' 

\x\\x\\x\\x\\ 



■ y y y- -r "^"V""^''-^-^ yyyyyyyyy 



Q^ {Gey/cf 



10 



Figure 14. Comparison of the twist analyses of the proton polarized Mn {Q"^) (a,c,e) and un- 
polarized M„ (Q^) (b,d,f) Nachtmann moments adopting the SGR technique resulting from this 
work and from Ref. |6[, respectively. The meaning of the lines and shaded areas is the same as in 
Fig. 13. Note that the scale of the vertical axis is different for polarized and unpolarized moments. 
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